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ABSTRACT 



We report the discovery, with NASA's Rossi X-ray Timing Explorer (RXTE), of the 
first sub-millisecond oscillation found in a celestial X-ray source. The quasi-periodic 
oscillations (QPO) come from Sco X-l and have a frequency of approximately 1100 Hz, 
amplitudes of 0.6-1.2% (rms) and are relatively coherent, with Q up to ~10 2 . The 
frequency of the QPO increases with accretion rate, rising from 1050 to 1130 Hz when 
the source moves from top to bottom along the normal branch in the X-ray color-color 
diagram, and shows a strong, approximately linear correlation with the frequency of the 
well-known 6-20 Hz normal/flaring branch QPO. We also report the discovery of QPO 
with a frequency near 800 Hz that occurs, simultaneously with the 1100 Hz QPO, in 
the upper normal branch. We discuss several possible interpretations, one involving a 
millisecond X-ray pulsar whose pulses we see reflected off accretion-flow inhomogeneities. 
Finally, we report the discovery of ~45Hz QPO, most prominent in the middle of the 
normal branch, which might be magnetospheric beat-frequency QPO. 

Subject headings: stars: individual (Sco X-l) — stars: neutron — pulsars: general 
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1. Introduction 

The characteristic time scale for motion of matter 
near a gravitating object, the dynamical time scale 
Tdyn, is (r 3 1 'GM) 1/>2 , where r is the distance to the 
object and M its mass. For r = R, the size of the 
object, Tdyn is relevant to possible spin and vibra- 
tion periods. For a neutron star Td yn extends down 
to values of well below 1 ms. Sub-millisecond vari- 
ability could therefore be produced in an accreting 
neutron star by various different mechanisms. In this 
paper, we present the first conclusive evidence for sub- 
millisecond variability in an accreting neutron star, 
ScoX-I. A preliminary announcement of this work 
was already made in Van der Klis et al. (1996). 

ScoX-1 (Giacconi et al. 1962), the brightest persis- 
tent X-ray source in the sky, is a Z source (Hasinger 
and Van der Klis 1989), a luminous low-mass X-ray 
binary containing a low-magnetic-field neutron star. 
In the X-ray color-color diagram it shows the lower 
two of the three branches in the canonical "Z track" , 
the normal branch (NB) and the flaring branch (FB) 
(Middlcditch and Priedhorsky 1986, Priedhorsky et 
al. 1986, Hertz et al. 1992, Dieters and Van der Klis 
1996). Its power spectrum is characterized by the 
presence of 1-5% amplitude, 6-20 Hz normal/flaring 
branch QPO (N/FBO). According to the standard in- 
terpretation (e.g., Hasinger et al. 1990, Lamb 1991), 
when the mass transfer rate M to the neutron star in- 
creases, the source moves down along the NB, passes 
through the "vertex" of the two branches, where it 
reaches the Eddington critical rate, and then moves 
up the FB. At the same time the QPO frequency in- 
creases from ~6 Hz on the NB to ^20 Hz on the lower 
FB. Further up the FB the QPO disappear. Below, 
we report three new, previously unknown QPO phe- 
nomena in this source, and describe their dependence 
on M. 

2. Observations 

ScoX-1 was observed three times with the propor- 
tional counter array (PCA) onboard NASA's Rossi X- 
ray Timing Explorer (RXTE; Bradt, Rothschild and 
Swank 1993), on 1996 Feb. 14 from 9:14 to 13:25 UT, 
Feb. 18 4:46-8:40 UT and Feb. 19 10:09-14:56 UT, here- 
after observations 1, 2 and 3, respectively. Each ob- 
servation covered three satellite orbits with ~2500 s of 
data, separated by intervals of ~1500s due to Earth 
occultation and/or South-Atlantic Anomaly passage. 
The count rate was ~ 10 5 c/s in observations 1 and 2. 



At the start of orbit 1 of observation 3 the count rate 
increased to ~1.45 10 5 c/s, which tripped the PCA 
high- voltage safety switch-off. RXTE was then moved 
slightly off ScoX-1 to ^50% collimator efficiency; 
later count rates were between 0.5 and 0.8 10 5 c/s. 
The <20keV background was ^90 c/s and has been 
negelected. 
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Fig. 1. — X-ray color-color diagram. Soft color is 
the (3-5)/(l-3) keV, hard color the (7-20)/(5-7) keV 
count-rate ratio. Each data point corresponds to 16 s 
of data; filled circles, open circles and crosses corre- 
spond to observations 1, 2 and 3, respectively. The 
statistical uncertainties are ^0.003 in soft color and 
~0.005 in hard color. The systematic uncertainties in 
Z track location between observations (see text) are 
-1%. 

During orbits 1 and 3 of each observation, data 
were collected with a time resolution of initially 0.5 ms 
and later, after discovery of the 1100 Hz QPO in ob- 
servation 1, usually 0.25 ms. During orbit 2 of each 
observation, the time resolution was set to 16 /US, and 
"double events" were additionally recorded at 64 /xs 
time resolution. Double events, two events detected 
within 6 /zs at two different anodes of a PCA detec- 
tor, are normally mostly due to charged particles. 
However, for the very high count rates from ScoX-1 
they are mostly due to two source photons. In the 
data analysis the double events, counted as two, were 
added to the single-photon data. This led to a con- 
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sidcrable increase in sensitivity to time variability. 
All high-timc-rcsolution data were recorded in the 2- 
20keV band. During the entire run 16-s resolution 
data were additionally recorded in 129 spectral chan- 
nels covering the 2-60 keV band. 

3. Analysis and Results 




500 1000 1500 2000 

FREQUENCY (Hz) 



Fig. 2. — Power spectra of combined single- and 
double-event data (see text) showing simultaneous 
800 and 1100 Hz QPO (top; orbit 2 of observation 
1), and a narrow 1100 Hz QPO peak (bottom; orbit 2 
of observation 2). The spectra are Leahy- normalized; 
the offsets of the continua from a level of 2.0 and the 
slope of the continuum in observation 1 are due to 
instrumental deadtime effects. 

We used the 16-s data to construct X-ray color- 
color diagrams. There was a considerable offset 
between the Z track seen in observations 1 and 2 
and that in observation 3, which we attribute to 



energy-dependent detector effects (collimator reflec- 
tions and/or rate-dependent gains). After an empiri- 
cal correction that assumes the Z track did not move 
between our observations (e.g., Dieters and Van der 
Klis 1996), a typical ScoX-1 color-color diagram re- 
sults (Fig. 1). During observation 1 the source moved 
gradually down the NB, during observation 2 it moved 
up and down between the lower NB and FB, and in 
observation 3 it moved gradually down the FB and 
into the vertex region. We assigned numbers Sz to 
positions on the Z track according to curve length 
along the track, where we set the top of our NB to 1, 
and the vertex to 2. 

Approximately 1100-Hz QPO (Fig. 2) were ob- 
served on all occasions that the source was on the 
NB or very low on the FB, and the observational set- 
up was adequate (full collimator effciency to ensure 
sufficient sensitivity, and sufficient time resolution). 
This was the case in orbit 2 of observation 1 and in 
orbits 1 and 2 of observation 2. The QPO were seen 
independently in single and in double-event data. Ta- 
ble 1 lists the results of functional fits of a constant 
plus a power law plus one or two Lorentzian peaks to 
power spectra of the high time resolution data. All 
fits were statistically acceptable with reduced \ 2 <~ 1. 
The QPO frequency, fnoo 5 increased with M, from 
1050 Hz at S z = 1.25 to 1130 Hz at S z = 2.1. Note 
that for lack of time resolution we have no information 
about the 1100 Hz QPO in the S z range 1.4-1.8. The 
QPO were sometimes quite coherent, with Q = vj Av 
up to 10 2 . Fractional rms amplitudes, corrected for 
differential deadtime (Van der Klis 1989) assuming a 
paralyzable process (Zhang et al. 1995) with a dead- 
time of 10 /xs (Zhang 1995) ranged between 0.9 and 
1.2% in observation 1 and 0.6-0.9% in observation 
2. There was no clear dependence on M within each 
observation. 

QPO with a frequency of approximately 800 Hz 
were observed simultaneously with the 1100-Hz QPO 
(Fig. 2 ( top )) only when the source was in the upper 
NB (Sz<1.36). From studying the short-term power- 
spectral variations we can exclude that the two peaks 
are due to one peak moving in frequency, unless it 
moves on time scales shorter than 32 s. In orbit 2 of 
observation 1 the QPO frequency, i^oo; generally in- 
creased (from ^800 to ^830 Hz) when Sz rose from 
1.25 to 1.35 (z^iioo increased from 1050 to 1075 Hz over 
the same range). The 800 Hz QPO were relatively 
broad, 50 to > 100 Hz and had amplitudes between 
0.9 and 1.2% (rms), decreasing slightly with Sz- In 
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orbit 1 of observation 1 (Sz<1.16) the sensitivity was 
less due to the lack of double-event data, but there 
was some evidence (3a) for the presence of a broad 
peak near 720 Hz. 

A third new QPO feature, located near 45 Hz, was 
detected on the NB. These QPO had a FWHM of 
typically 15-25 Hz, an amplitude of ~1% (rms) and 
were most prominent on the middle of the NB, near 
Sz = 1.5. They were detected all the way from Sz 
= 1.1 to 1.9 with little change in frequency. Fig. 3 
shows the power spectrum of data obtained in orbit 
3 of observation 1, with 6 Hz QPO (the N/FBO) and 
45 Hz QPO. 
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Fig. 3. — Power spectrum of 0.5 ms resolution data 
obtained in orbit 3 of observation 1, showing a strong 
N/FBO peak centered on 6 Hz, and a weaker feature 
near 45 Hz. The spectrum is Leahy-normalized; the 
offset of the continuum from a level of 2.0 is due to 
dcadtimc effects. 



The dependencies of i^noo and the N/FBO fre- 
quency ^ N/FBO on Sz are very similar. The relation 
between vuoo and ^ N/FBO is plotted in Fig. 4. There is 
a strong correlation. The relation fits a straight line, 
with fnoo = 7.3 
well. 
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Fig. 4. — Relation between the frequencies of the 
1100 Hz QPO and the N/FBO as derived from data 
where these two phenomena were detected simultane- 
ously. A strong correlation is evident. 



4. Discussion 

The discovery of near-millisecond (Strohmaycr, Zhang 
and Swank 1996, Van Paradijs et al. 1996, this pa- 
per) and submillisecond (Van der Klis et al. 1996, this 
paper) oscillations from accreting low-magnetic-field 
neutron stars was expected on the basis of general ar- 
guments (see Introduction). To find and study such 
phenomena was one of the primary objectives of the 
RXTE mission. There is a range of phenomena that 
could take place on these time scales, and we shall not 
attempt to be exhaustive here. The fact that both 
^noo and vsoo increase considerably as a function of 
M suggests an origin in accretion phenomena rather 
than neutron star vibrations (e.g., McDermott, Van 
Horn and Hansen, 1988). Therefore, we concentrate 
on that class of models. 

Disk frequency. The Kepler frequency vk at the 
inner edge of the accretion disk, bounded by a small 
(16-20 km) magnetosphere is in the correct range to 
be identified with either ^noo or fsoo- It would in- 
crease with M, as observed. The M dependence of 
the 800 Hz QPO, strongest at low M and disappear- 
ing towards higher M is reminiscent of the "horizontal 
branch QPO" (HBO) in other Z sources (Van der Klis 
et al. 1985, Van der Klis 1995 for a review). These 
HBO, in the beat-frequency model (Alpar and Sha- 
ham 1985, Lamb et al. 1985), are caused by "clumps" 
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in the inner disk, which are most prominent at the 
lowest M levels. The 800 Hz QPO (or, alternatively, 
the 1100 Hz QPO) could be a direct signal from these 
clumps, either by reflection or obscuration of the X- 
ray flux from the central regions, or by interaction 
with azimuthal structure in the radial flow (see be- 
low). 

Beat frequency. As no HBO have been detected 
as yet from ScoX-1, an alternative possibility is that 
the 800 Hz QPO are HBO, and their frequency the 
beat frequency vb = u{vk — v s ) between vk and 
the neutron star spin frequency v s , where n is the 
field symmetry factor. This would make the HBO in 
ScoX-1 much faster than those in other Z sources, 
whose frequency is 15-55 Hz, and would require a rel- 
atively slow neutron star spin rate (not more than a 
few 100 Hz). When this Letter was about to be sub- 
mitted, the detection was reported in 4U 1728—34 of 
363 Hz oscillations during X-ray bursts (Strohmayer 
et al. 1996) which were interpreted as due to the 
neutron star spin, and two variable-frequency QPO 
peaks (one of which the ^800 Hz QPO reported by 
Strohmayer, Zhang and Swank, 1996) whose separa- 
tion remained near 363 Hz. Clearly, this suggests a 
beat frequency model interpretation wit vb, vk an d 
v s all visible. In the case of ScoX-1, the difference 
^noo — ^800 is consistent with being constant near 
250 Hz, but the range in frequencies over which both 
peaks are detected simultaneously is too small to draw 
any strong conclusion from this; there is no evidence 
for an oscillation near 250 Hz. 

The new 45 Hz QPO in ScoX-1 are in the range 
of HBO frequencies seen in the other Z sources. The 
dependence of their amplitude on Sz is clearly differ- 
ent from that of HBO in GX5-1 and CygX-2, which 
are strongest at the top of the NB (and in the HB), 
not in the middle of the NB. In GX 17+2, 60 Hz QPO 
have been discovered in the NB which may be a sim- 
ilar phenomenon (Wijnands ct al. 1996, in prep.). If 
the 45 Hz QPO in ScoX-1 are beat-frequency QPO, 
then from the fact that its frequency is approximately 
constant we conclude that the 1100 Hz and the 800 Hz 
QPO, whose frequencies vary by tens of Hz, can not 
be the disk frequency 

Doppler- shifted pulsar. The strong correlation we 
observe between ^noo and f N/FBO is not explained by 
any of the models discussed above. It could in princi- 
ple arise because both frequencies vary as a function 
of M independently of each other, but this seems un- 
likely. The ^ N/FBO vs. Sz relation is known to jump at 



Sz=2 (see, e.g., Van der Klis 1995), and ^noo should 
then, independently, show a similar jump in its de- 
pendence on Sz to maintain the correlation. It seems 
more likely that one frequency is by some mechanism 
directly derived from the other one. 

We briefly explore a mechanism where this is the 
case. In the Fortner, Lamb and Miller (1989) model 
for N/FBO, part of the accretion takes place by way 
of an approximately radial inflow, and blobs (inho- 
mogeneities) in this flow produce the N/FBO. For 
a constant-size radial flow, ^ N/FBO is proportional to 
the inflow velocity v r . If the 1100 Hz QPO are the 
Dopplcr-shiftcd pulsar signal, which we see reflected 
off blobs in the flow, then when v r varies, vuoo would 
vary linearly with ^ N/FBO , as observed. For a blob 
moving towards the star along a radius vector that 
makes an angle <j> with the line of sight, ^noo would 
be proportional to (1 + (1 — cos(f>)v r /c)v pu i s , where 
v V uis is the pulse frequency. Radial flow velocities 
up to at least 0.04c, and spin frequencies of ~500 or 
~1000Hz, depending on whether one or two magnetic 
poles contribute to the signal, are required for this 
model to work. To get the relatively coherent QPO 
signal we observe, the range of angles 4> contributing 
to the signal should be small. This would be true, 
for example, if the magnetic and rotation axes are 
nearly aligned. Aligned rotation would also explain 
why millisecond pulsars in low-mass X-ray binaries 
have been hard to find. This model would require a 
spin period ~lms, faster than any neutron star spin 
yet measured, as only one pole would contribute (the 
other pole could produce another QPO peak). 

Note. Further RXTE observations of ScoX-1 on 
May 24-25 again show the two high-frequency QPO 
peaks, this time with frequencies near 900 and 600 Hz. 
Their separation is significantly (>10cr) larger than in 
the February observations. This excludes the beat- 
frequency interpretation where the peak separation 
is predicted to be constant at the neutron star spin 
frequency. 
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Table 1 
qpo parameters 



1100 Hz 



800 Hz 



N/FBO 



Frequency FWHM 
(Hz) (Hz) 



Frequency 
(Hz) 



FWHM 

(Hz) 



Frequency 
(Hz) 



FWHM 

(Hz) 



Observation 1 orbit 2. 



1.252±0.027 


1050.3±1.1 


41±5 


803±3 


88±11 


4.9±0. 


6 


13±3 


1.282±0.040 


1069.8±8.2 


104±30 


794±20 


145 +189 
140 -50 


6.06±0. 


25 


9.3±1.1 


1.336±0.032 


1080.0±4.8 


50±19 


833±12 


80 fix 


5.74±0. 


28 


10.7±1.1 


1.340±0.050 


1072.5±5.4 


78±18 


829±5 


49±19 


6.06±0. 


22 


8.5±0.8 


1.357±0.022 


1069±15 


178+^f 


820±12 


80 fix 


5.64±0. 


27 


9.9±1.1 



Observation 2 orbit 2. 
1.866±0.022 1091. 4±1. 8 
1.929±0.050 1100.8±2.0 

Observation 2 orbit 1. 
1.966±0.032 1090.1±2.3 



2.057±0.006 
2.092±0.025 



1129.7±2.3 
1135±13 



22±5 
18±4 



18±6 

a +88 
-40 



68^ 



7.85±0.04 
9.42±0.07 



7.66±0.09 
13.64±0.17 
13.46±0.04 



4.44±0.12 
7.24±0.23 



5.81±0.27 

6.6±0.6 
9.60±0.18 



Fit function: sec text. All errors correspond to unreduced A\ 2 = 1. Fit range, frequency resolution 
and number of degrees of freedom for the 1100 and 800 Hz QPO fits were 200-2000 Hz, 10 Hz and 171 
for observation 1 orbit 2, 512-1536 Hz, 5 Hz and 200 for observation 2 orbit 1 and 200-2000 Hz, 5 Hz 
and 356 for observation 2 orbit 2, respectively. 
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